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ABSTRACT 

We explore the formation of secondary Population III (Pop III) stars under radiation hydrodynamic 
(RHD) feedback by a preformed massive star. To properly treat RHD feedback, we perform three- 
dimensional RHD simulations incorporating the radiative transfer of ionizing photons as well as H2 
dissociating photons from a preformed star. A collapsing gas cloud is settled at a given distance from 
a 120M Q Pop III star, and the evolution of the cloud is pursued including RHD feedback. We derive 
the threshold density depending on the distance, above which the cloud can keep collapsing owing 
to the shielding of H2 dissociating radiation. We find that an H2 shell formed ahead of an ionizing 
front works effectively to shield the H2 dissociating radiation, leading to the positive feedback for 
the secondary Pop III star formation. Also, near the threshold density, the envelope of gas cloud is 
stripped significantly by a shock associated with an ionizing front. By comparing the mass accretion 
timescale with the Kelvin-Helmholtz timescale, we estimate the mass of secondary Pop III stars. It 
turns out that the stripping by a shock can reduce the mass of secondary Pop III stars down to 
w 20M Q . 

Subject headings: theory:early universe — galaxies: formation — radiative transfer — molecular 
processes — hydrodynamics 



1. INTRODUCTION 

The formation of Population III (hereafter 
Pop III) stars has bee n explored extensively 

by many authors (iBromm, Coppi fc Larson] 

1999L 12001 iNakamura fc Umemur al 119991 l27)0lf 
Abel. Brvan. Normanl I2000L [200l lYoshidal [20061 : 
Yoshida et alj 120061: lO'Shea fc Norman 200fl7 and it 



is expected that the mass function of Pop III stars 
is top-heavy with the peak of 100 — 1000M©. But, 
the mass range of Pop III stars is not so clear. The 
enhanced H2 cooling in pre - ionize d gas may r educe 
the mass (IShapiro fe KangJ 11987k ISusa et all 119981 : 
lOh fc Haiman I [2002T l . If the HD cooling is efficient in 
fossil HII regions, the mass could be some 10M P , (e.g. , 
Uehara & Inutsuka 2000; Nakamura & Umemura 2002; 
Naga kura fc Omukail 120051: iJohnson fc Bromm 120061: 
Greif fc Bromml 120061 : lYoshida. Omukai fc HernquistJ 
2007[) . Also, the variations of cosmological density 



fluctuations may lead to the format ion of less massive 
Pop III stars with some 10M S (lO'Shea fc Normanl 
2007). However, the radiation hydrodynamic effect on 
the mass of secondary Pop III stars forming nearby 
a preformed massive star has not hitherto explored 
in detail. From an observational point of view, the 
elemental abundance patterns of hyper-metal-poor stars 
well match the yields by supernova explosions with a 
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progenitor mass of ~ 25M Q ([Umeda fc Nomotol 120031 : 
llwamoto et al.|[2005f ). This implies that Pop III stars of 
some 10M Q are likely to have formed before the first 
metal enrichment. 

Metal-free first objects are expected to collapse at 
20 < z < 30, forming a minihalo wit h the mass of 
» lC ^ Mm and the e xtent o f » lOOpc (iTegmark et all 
19971 iNishi fc Susal H99l iFuller fc Couchmanl |2005T 
Yoshida et al. I l2003f h In the course of bottom- up 
structure formation, such minihaloes merge to form first 
galaxies at z > 10, having virial temperature > 10 4 K 
and mass > 1O 8 M . Even in the evolution of first 
galaxies, Pop III stars can play a significant role, since 
an appreciable number of stars c an form from metal-free 
component in inte rstellar gas (|Tornatore et all 120071 : 
IJohnson et al|[2b08f ). Also, those Pop III stars are l ikely 
to be r esponsible for the reioniza t ion of t he universe (ICen 
2003) ICiardi Ferrara. fc White! [200l IWvithe fc Loet 
2004 ISomerville fc Livid l2003t TSokasian et alj 12004 
Murakami et all [2005), and the metal enrichment 
of interg alactic medium (INakamura fc Umemural 
20011; IScannapieco. Ferrara. fc Madaul 120021 : 



Ricotti fc OstrikerT 12004) thro u gh supernova (SN) 
dMori et aTT 



KiL_ . 

explosions dMori et all 20021: IBromm et all 120031 : 
iKitavama fc Yoshidal 120051 : iGreif et all l2007f l Hence, 
the mass range of Pop III stars affects the galaxy 
formation and the evolution of intergalactic matter. 

If first stars are very massive, they emit intensive 
ultraviolet (UV) radiation, which significantly influ- 
ences the subsequent star formation in minihaloes or 



first g alaxies (lHaiman et al.l 1997; Omukai & Nishi 
19991: IHaiman et al l 120001 blover fc Brandl 120011 : 



Machacek et al.ll2001l ). In a Pop III minihalo, primordial 



fluctuations generate density peaks owing to gravita- 
tional instability induced by H2 cooling. A highest peak 
col lapses earlier to form a first Pop III star as shown 
by lAbel. Bryan, Normanl ((2002). Subsequently, lower 
peaks collapse to form cloud cores with the density 
of 10 2 - 4 cm- 3 and the temp erature of 200 - 500K 
(|Bromm. Coppi fc Larsonll2002f ) . Hence, later collapsing 
cores could be affected by the radiative feedback by the 
first star. In first galaxies, the interstellar medium can 
be pre-ionized by shock, since the virial temperature is 
higher than 10 4 K. The pre- ionization c an promote rapid 
H2 fo r mation (|Shapiro fc Kang 1119871 iKang fc Shapiro I 
119921 iSusa et alJ 119981 : lOh fc Haiman I |200^ T As a 
result, the temperature of interstellar gas in first 
galaxies can be lowered down to 100K through efficient 
H 2 cooling (|Johnson et al.l l2008h . Low temperature 
collapsing clouds in first galaxies could be subject to the 
UV feedback by Pop III stars. 

If a first star is distant by more than lpc, dense 
clouds are r eadily self-shielded f r om th e ultraviolet (UV) 
radia ti on dTaiiri fc Umemural 119981 : iKitavama et al 
2001b ISusa fc Umemural l2004allbl; IKitavama et al. 
2004: iDijkstra et alJ 12004 lAlvarez et all 12009) . Thus, 
the photoevaporation by UV heating is unlikely to 
work devastatingly. However, Lyman- Werner (LW) 
band radiation (11. 18-13. 6eV) photodissociates H2 
molecules, and therefore can preclude a cloud core from 
collapsing. Hence, the photodissociatio n of H2 may 
lead t o momentous neg a tive f ee dback (IHaiman et al.l 
19971 lOnTukai fc NishiT [1991 [Haiman et all 120001 : 
Glover fc Brandl 120011 : iMachacek et al.ll2001h . Recently, 
Susa I (|2007l ) performed 3D radiation hydrodynamic 
simulations on the photodissociation feedback, and 
derived the criteria of the feedback. 

However, ionizing radiation (>13.6eV) is related to 
H2 formation in an intricate fashion. An ionization 
front (I-front) driven by ionizing radiation propagates in 
a collapsing core. The enhanc ed fraction of electrons 
catalyzes rapid H 2 fo r mation ([Shapiro fc Kang I [19871: 
Kane fc Shapiro 1119921: ISusa et al.lll998l : lOh fc Haiman I 
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In particular, the mild ionization ahead of the I- 
front can generate an H 2 shell, which potentially shield s 
H2 dissociating photons (|Ricotti. Gnedin. fc S hull 2001). 
This mechanism is likely to work positively to form Pop 
III stars. When UV irradiates a dense core, the I-front 
changes from R-type on the surface to D-type inside the 
core. The transition occurs via an intermediate type 
(M-typ e), which is accompanied with the generation of 
shock (|Kahnlll954l ). The shock can affect significantly 
the collapse of the core. This is a totally radiation hydro- 
dynamic (RHD) process. Such radiation hydrodynamic 
feedback has been investigated by ID spherical RHD sim- 
ulations (lAhn fc Shapiroll2007D . 2D cylindrical RHD sim- 
ulations (IWhalen et al.l[2008h . and 3D RHD simulations 
(jSusa fc Umemurall2006| ). The results by 2D and 3D sim- 
ulations are in good agreement with each other. How- 
ever, quantitative feedback effects, including the com- 
pression and stripping by a shock, have not been explored 
yet over a wide parameter range in 3D simulations. 

In this paper, we attempt to investigate quantita- 
tively the RHD feedback on a collapsing cloud by a pre- 



formed first star. For the purpose, we solve the three- 
dimensional RHD in a wide range of parameter space. 
In the simulations, the radiative transfer of ionizing pho- 
tons as well as H2 dissociating photons is self-consistently 
coupled with SPH hydrodynamics. In §2, the numerical 
method is described. The setup of simulations is pre- 
sented in §3, and numerical results are given in §4. In 
§5, the formation criteria of secondary Pop III stars are 
numerically derived taking RHD feedback into account. 
In §6, the mass of secondary Pop III stars is estimated. 
§7 is devoted to the conclusions and discussion. 

2. NUMERICAL METHOD 

We perform numerical simulations wit h the radia tion 
hydrodynamics code that we developed ([Susa 200j|). In 
this section, we briefly summarize the scheme. The code 
is designed to investigate the formation and evolution of 
the first-generation objects at z > 10, where the radiative 
feedback from various sources plays important roles. The 
code can compute the fraction of the chemical species 
e~, H + , H, H~, H 2 , and Hj through fully implicit time 
integration. It also can deal with multiple sources of 
ionizing radiation, as well as the radiation at the Lyman- 
Werner band. Although there may be heavy elements 
in first galaxies through Pop III supernova explosions 
therein, the metallic cooling is not d ominant as long as 
the m etallicity is lower than 10~ 2 Z Q (jSusa fc Umemural 
2000). In this paper, we ignore the metallic cooling. 

The hydrodynamics is calculated by the smoothed par- 
ticle hydro dynamics ( S PH) method. We use the version 
of SPH by Umemura (1993 ), w i th the modification ac- 
cording to Isteinmetz fc Mullerl (I1993D. We adop t the 
particle resizing formalism bv lThacker "elaTI ((20001) . 

The non-equilibrium chemistry and radiative cooling 
for pri mordial gas are cal c ulated with the code devel- 
oped bv lSusa fc Kitavamal (|2000l ), where the H 2 cooling 
and re action rates are mostly taken from lGalli fc Palla 
(1998). As for the photoionization process, we em- 
ploy the so-called on-the-spot approximation (|Spitzerl 
1978), where it is assumed that recombination photons 
to ground states of hydrogen are absorbed promptly on 
the spot, while recombination photons to excited states 
can escape the medium. In the present version of our 
code, we do not take into account helium ionization. For 
pure hydrogen gas with on-the-spot approximation, we 
do not need the frequency bins for the radiation field. 
All we need is the optical depth at Lyman limit, since 
we know the frequency dependence of the optical depth. 
We can obtain the "precise" photoionization rate or pho- 
toheating rate as functions of the o ptical depth at Lyman 
limit before we start simulations (jSusa I [2006). We also 
remark tha t in case we include helium, we nee d more 
frequ encies (|Nakamoto. Umemura. fc Susa I 120011 : ISusa I 
[2009) . 

The optical depth is integrated using the neighbor lists 
of SPH particles. It is sim ilar to the code described in 
ISusa fc Umemura! (|2004al) . but now we can deal with 
multiple point sources. In our new scheme, we do not 
create as m any grid points on the li ght ray as the pre- 
vious code (jSusa fc Umemural l2004a[ ) does. Instead, we 
just create one grid point per SPH particle located in its 
neighbor. We find the "upstream" particle for each SPH 
particle on its line of sight to the source. Then the opti- 
cal depth from the source to the SPH particle is obtained 
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by summing up the optical depth at the "upstream" par- 
ticle and the differential optical depth between the two 
particles. 

The opacity against LW-band flux (Jlw) is calculated 
with the self-shielding function by iDraine fc Bertoldil 
(1996), which is given by 

where FLw.o is the incident flux, -ZVh 2 ,14 = 
■/VH 2 /10 14 cm -2 is the normalized H2 column den- 
sity, and 

/*(*)={ iv, x x i\ (2) 

The column density of H2 is evaluated by the method 
described above. We assume the same absorption for 
the photons within the LW-band between 11.26 and 13.6 
eV. The use of the self-shielding function ignores the ef- 
fects of Doppler shifts of LW lines. If the line is shifted 
more than thermal broadening width, the opacity at LW 
band is greatly reduced. In the present calculations, the 
collapsing gas cloud converges to Larson-Penston type 
similarity solutions. Thus, the infall velocity in the enve- 
lope of cloud is as large as the sound velocity. Moreover, 
the velocity around the core is much smaller than sound 
velocity. Therefore, we do not include this effect in the 
present simulations. 

The code is already parallelized with the MPI library. 
The computational domain is divided by the so-calle d 
orthogonal recursive bisection method (Dubinski 1996). 
The parallclization method for the radiation transfer part 
is similar to that of the multiple wave front me thod 
developed bvlNakamoto. Umemura. fc Susal (|2001[ ) and 
iHeinemann et a l. (2006), but it is changed to fit the SPH 
scheme. The co de is able to handle self-gravity with a 
Barnes-Hut tree (|Barnes fc HutllT986l ). which is also par- 
allcliz ed. The co de is tested for various standard prob- 
lems (|Susa 1120061 ) . We also take part in the code compar- 
ison project with other radiation hydrodynamics codes 
(llliev et all l200l [20091. and we find reasonable agree- 
ments with each other. 

The present simulations are mainly carried out with a 
novel hybrid computer system in University of Tsukuba, 
called FIRST simulator, which has been designed to sim- 
ulate multi-component self-gravitating radiation hydro- 
dyna mic systems with high accuracy (|Umemura et all 
I2007D . The FIRST simulator is composed of 256 nodes 
with dual Xeon processors, and each node possesses a 
Blade-GRAPE board, that is, the accelerator of gravity 
calculations. The peak performance is 36.1 Tflops. 

3. SETUP OF SIMULATIONS 

We use TV = 524, 288(= 2 19 ) SPH particles to perform 
radiation hydrodynamic simulations. We suppose a pri- 
mordial gas cloud undergoing run-away collapse nearby a 
preformed massive star. Dark matter is neglected in the 
present calculations, because we are interested in rather 
dense cores with n > 10 2-4 cm -3 that start to become 
baryon-dominated (e.g. Abel, Bryan & Norman 2002). 
The tidal field from other halos could have some effects 
on the dynamics of collapsing gas, which is not taken into 
account in the present calculation. The chemical com- 
positions are i nitially assumed to h ave the cosmological 
residual value (jGalh fc Palla II1998D . The mass of cloud 



is Mb = 8.3 x 10 4 Af Q in baryonic component. The gas of 
the cloud is initially at rest and have "top-hat" density 
distribution. We also assume initially uniform temper- 
ature. Remark that the density distribution converges 
to ordinary p cx r~ 2 profile as the collapse proceeds 
(e.g. Abel, Bryan & Norman 2002). Also, we distribute 
a low-density (n ~ 10 _1 cm~ 3 ) uniform gas around the 
collapsing cloud. When the central density of cloud ex- 
ceeds a certain value, n on , we ignite a 120M Q Pop III 
star. The luminosity and the effective temperature of the 
source star are taken from iBaraffc. He ger fc Wooselvl 
(200l|). The ionizing photon number per unit time is 
A^ion = 1-3 x 10 50 s _1 , whereas the LW photon number 
luminosity is A LW = 1.4 x 10 49 s _1 . The source star is 
located at the distance D from the center of collapsing 
cloud. 

First, we see fundamental physical processes of radia- 
tive feedback using high/low temperature core models. 
Then we derive the radiation hydrodynamic feedback 
criteria by changing n on and D. The basic difference 
between two models is the initial temperature of the 
clouds. The high initial temperature results in the low 
temperature core (LTC), while the low initial tempera- 
ture does in the high temperature core (HTC). In Fig. 
[U the cloud collapse in two cases is shown by tempera- 
ture versus density at the center. In a LTC model, the 
relatively high initial temperature (Tj n i = 350K) allows 
the efficient formation of H2 molecules and also leads 
to the relatively slow gravitational contraction. In this 
model, the ratio between gravitational energy to initial 
thermal energy |PF|/£7 ~ 2. Consequently, the core tem- 
perature goes down quickly to w 150K at the density 
Tin ~ 100cm~ 3 . Since the size of core is approximately 
the Jeans length, a smaller core forms in a LTC model. 
It is worth noting that the density distribution of collaps- 
ing gas for riH > 10 3 cm~ 3 is close to the one obtained 
in c osmological simulations including dar k matter grav- 
ity (|Q'Shea fc Normanl [2007t [Su^a1l200l . In an HTC 
model, the formation of H2 molecules is not so rapid ow- 
ing to the lower initial temperature (Ti n ; = 100K). Also, 
the initial low pressure (|W|/Z7 ~ 3 — 4) allows the rapid 
gravitational contraction. Hence, the temperature in a 
collapsing phase becomes higher, and eventually a larger 
core forms. The difference of core size is significant for 
the radiative feedback. 

The basic models investigated in this paper are shown 
in Table CD LTC-ION and HTC-ION include both ioniz- 
ing and LW-band radiation from a 120M Q star respec- 
tively for LTC and HTC models, while LTC-LW and 
HTC-LW include only LW-band radiation respectively 
for LTC and HTC models. T c m i n is the minimum core 
temperature after the H2 cooling instability. In these 
models, D = 40pc and n OD = 10 3 cm~ 3 are assumed, but 
D and changed over a wide range to derive the 

feedback criteria. 

4. NUMERICAL RESULTS 

Here, we see the fundamental physical processes for 
four basic models shown in Table [TJ 

4.1. Evolution of Cloud Cores 

The time evolutions of cloud cores for four models are 
shown in Fig. [5] In upper two panels the tempera- 
ture and H2 molecule abundance at the cores are shown 
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Fig. 1. — Cloud collapse without radiative feedback. The tem- 
perature at the center is shown against the central density for the 
initial temperature of 350K (red curve) and 100K (blue curve). 



TABLE 1 
Basic Models 



2Myr (the second peak in a red solid curve in the up- 
per left panel in Fig[2]) raises the temperature quickly, 
so that the core cannot be gravitationally bounded and 
bounces eventually. Such intricate behavior of a shock is 
explained further in the following. 

4.2. Convergence check 

Ideally, it is appropriate to perform convergence check 
for all the models, however, we check for HTC-ION 
model only, because of the limited computational time. 
We run three additional simulations with particle number 
N = 2 17 , 2 18 , 2 20 , while the canonical run uses N = 2 19 
particles. Since the highest density peak collapses and 
we are now interested in the collapse criteria, we plot 
the evolution of the density peak on density-temperature 
plane and density-H2 abundance (j/h 2 ) plane for four dif- 
ferent resolutions. As shown in the Figure [3l all of the 
results agrees fairly well, although the lowest resolution 
case (long dashed line) shows some deviation from oth- 
ers. Thus, the the resolution of the present simulations 
are enough to describe the collapse criteria. 

4.3. Cloud Structure 

The evolution shown above indicates that ionizing ra- 
diation works to promote H2 formation, and therefore 
help the cloud to collapse. In Fig. [4] the spatial dis- 
tributions of H2 molecules on the plane including the 
.sym metry axis are shown for LTC-ION and HTC-ION 



Model 


T c , min [K] 


Ionizing Flux 


LW-band 


D [pc] 


LTC-ION 


150 


Yes 


Yes 


40 


LTC-LW 


150 


No 
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40 


HTC-ION 


300 


Yes 


Yes 


40 


HTC-LW 


300 


No 


Yes 
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[cmirljjd cls. Two snapshots are taken at the epoch of lMyr 
and 2Myr after UV irradiation. FigsE] and [6] show the 
gas number density and temperature distributions. Also, 
we plot various physical quantities along the symmetry 
■a*is for all four models in Fig [7] 
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against the peak density for LTC-ION and HTC-ION 
models, while in lower two panels they are shown for 
LTC-LW and HTC-LW models. 

The cloud evolution is regulated by the formation of 
H2 molecule, which is an indispensable coolant for cloud 
collapse. As seen in Fig. [2j the H2 abundance drops 
down to a level of 10~ 9 — 10~ 6 immediately by UV ra- 
diation after the ignition of a source star (right panels, 
where the ignition time is marked as OMyr). Then, owing 
to the inefficiency of H2 cooling, the gravitational con- 
traction becomes adiabatic and temperature rises until 
lMyr (left panels). During the adiabatic contraction, H2 
molecules retrieve, since the H2 formation rate is raised 
in adiabatically heated cores. 

In HTC-ION model, the heating by a shock associated 
with an ionization front helps to increase the tempera- 
ture. As a result, the H2 formation is restored and H2 
fraction exceeds ~ 10~ 4 . Thus, the core can cool and 
keep collapsing. In LTC-LW, and HTC-LW models, the 
H2 abundance drops down to a level lower than 10~ 8 
after UV irradiation, and the temperature continues to 
increase owing to insufficient cooling. Although H2 for- 
mation is recovered, H2 fractions are not sufficient for 
the cloud to keep collapsing. Eventually, the cloud core 
bounces. In LTC-ION model, the physical processes are 
more complicated. The shock heating before lMyr works 
importantly to restore H 2 formation. Resultantly, the H 2 
abundance reaches to a level of ~ 10~ 5 — 1CV 4 . However, 
the second hit by a diffracted shock between lMyr and 



In LTC-ION model (upper two panels in Figs[4][6]), 
the core is self-shielded against ionizing photons and a 
bow shock forms preceding an ionization front. The 
bow shock strips the envelope of a collapsing cloud. H2 
molecules are produced efficiently ahead of an ionization 
front owing to the enhanced catalytic effect of free elec- 
trons (see also Fig[7]). As a result, an H2 shell is built 
up and it shields effectively the H2 dissociating radiation 
from a source star. In Fig. [5J H2 column density from 
the source star is plotted at lMyr after UV irradiation. 
In fact, Nn 2 at the center of cloud is several hundred 
times higher than the critical value 10 14 cm~ 2 for the 
shielding of LW-band radiation. Thus, H2 dissociating 
radiation is significantly absorbed before it reaches the 
cloud center. However, as shown in the top-right panel 
in FigsHJini a diffracted shock front hits the cloud core 
before undergoing the run-away collapse. Consequently, 
the core bounces owing to the thermal pressure enhanced 
by the shock. It is also worth noting that we find "an H2- 
enriched tail" behind the core. This structure is formed 
by the collision of shock which is diffracted around the 
dense core. The temperature of the shock heated gas is 
high enough to feed electrons that catalyze the H2 forma- 
tion. (This structure could be interesting if it fragments 
into small clumps and yields protostars, although we do 
not find any sign of such a process in the present simula- 
tions. Regarding this issue, higher resolution simulations 
may be required.) 

In HTC-ION model (lower two panels in Figs[4][6]), a 
bow shock forms, similar to LTC-ION model. H2 column 
density at an H2 shell is comparable to that in LTC-ION 
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Fig. 2. — The cloud evolution with radiative feedback for the models listed in Table [T] The time variations of temperature (left two 
panels) and H2 molecule abundance (right two) at the cores are shown against the peak density. A red solid line depicts LTC-ION model, 
a blue solid line does HTC-ION (upper two panels), a red dotted line does LTC-LW, and a blue dotted line does HTC-LW (lower two 
panels). On each curve three stages are indicated by dots, which correspond to Oyr, lMyr, and 2Myrs after the ignition of a source star. 



model. However, in this case, the core size is larger than 
LTC-ION model. Thus, H2 column density at the cloud 
center is roughly ten times higher than that in LTC- 
ION model (see Fig. [8]). Also, a diffracted shock is not 
converged into the cloud core before the central density 
peak goes into run-away collapse phase. As a result, 
the cloud core can keep collapsing, although the cloud 
envelope is stripped by a bow shock. 

In LTC-LW model, a smaller core forms owing to its 
lower temperature. Since the column density of H2 
molecules inside the cloud is low, the self-shielding of 
H2 dissociating radiation is very weak (see Fig[7|). As 
a result, the photodissociation of H2 molecules is so in- 
tense to reduce the H 2 abundance to j/H2 ~ 1CP 8 . Con- 
sequently, the cloud core cannot cool and collapse. 

In HTC-LW model, although a large core forms owing 
to its higher temperature, the dissociating radiation from 
a source star destructs H2 molecules to a level of H2 
abundance as ?/H2 ~ 1CP 6 . Hence, the H2 cooling is not 
sufficient and the cloud core cannot collapse. 

To conclude, within these four models, only HTC-ION 
model succeeds in collapsing. The other models result 
in bouncing after UV irradiation. In other words, the 
positive feedback overwhelms the negative ones in HTC- 



ION model, whereas the negative feedback hinders cloud 
collapse in LTC-ION, LTC-LW, and HTC-LW models. 
However, the fate of clouds depends on the distance D 
from a source star and the core density n on at UV irradi- 
ation. Hence, we investigate a wider parameter space to 
derive the formation criteria of secondary Pop III stars. 

5. FORMATION CRITERIA 

In addition to four runs described in the previous sec- 
tion, we perform numerical simulations with various D 
and n on . Figs. [91 and [TTJ1 show the results of various runs 
for LTC and HTC models, respectively. The horizontal 
axis represents the density n on at turning on a source star 
and the vertical axis shows the distance D between the 
source star and the cloud center. In the runs denoted by 
crosses, the clouds cannot collapse until the end of simu- 
lations, i.e. 10 7 yr after the ignition of a source star. The 
open circles represent the runs with successful collapse. 
These runs can collapse even without ionizing radiation. 
On the other hand, the runs denoted by open triangles 
can collapse only if ionizing radiation is present. Here, "a 
collapsed cloud" is defined as the one whose central den- 
sity exceeds the numerical resolution limit determined by 
the Jeans condition. The density corresponding to the 
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Fig. 3. — Convergence check on density-temperature plane (left panel) and density-H2 abundance (yu 2 ) plane for HTC-ION model. In 
both of the panels, results from four runs are plotted. Short dashed: N = 2 20 , thick solid: TV = 2 19 (canonical), dashed: N = 2 18 , long 
dashed: N = 2 17 . 



resolution limit is ~ 10 6 — 10 7 cm~ 3 in the present sim- 
ulations. Solid curves show the collapse criterion in the 
presence of H 2 dissociating r adiation only , which is the 
same as the one obtained by iSusa I (|2007f ). Figs. [9] and 
[10] give the threshold density depending on the distance, 
above which the cloud can keep collapsing owing to the 
shielding of H2 dissociating radiation. 

Interestingly, there is no case where a cloud that can 
collapse without ionizing radiation fails to collapse by 
adding ionizing radiation. Near the threshold density, 
ionizing radiation can rather assist the collapse of a cloud 
that cannot collapse by H2 dissociating radiation, as 
shown by open triangles in Figs. [5J and [JO] Thus, we 
can conclude that ionizing flux alleviates the negative 
feedback by photodissociating radiation to some extent. 
In LTC models, the parameter region of open triangles 
is narrower, compared to HTC models. The contrast ba- 
sically comes from the difference in core temperature. 
In both models, the H2 shell shielding is almost the 
same, as shown in FigJH] Also, the mild shock preced- 
ing the ionization front heats up the core, which leads 
to efficient H2 formation and radiative cooling. This ef- 
fect is also reported by one-dimentional RHD simulations 
(]Ahn fe Shapiro 2007). However, the hydrodynamic evo- 
lution in the two models is very different. In HTC mod- 
els, the central part of the core can cool and collapse 
with aid of the mild shock heating as well as the H2 shell 
shielding. In LTC models, since the core mass and radius 
are smaller, a diffracted shock front is converged into the 
core (upper right panels in Figs[4][6]) before the collapse 
is completed. As a result, the central part of the core is 
heated up again. This second shock heating is too strong 
to keep the gas cloud bounded by self-gravity. Thus, the 
core bounces and evaporates after this second impact by 
a shock. This negative shock heating mechanism cancels 
the positive effects by the H2 shell shielding and mild 
shock heating in LTC-ION models. 

6. MASS OF SECONDARY POP III STARS 

To estimate the final stellar mass after the radiation 
hydrodynamic feedback, the mass accretion timescale 



is compared to the Kelvin-Hclmho ltz (KH) contrac- 
tion ti mescale, which is adopted from lO'Shea fc Normanl 

Fig. [UJ shows the mass accretion timescale as a func- 
tion of enclosed mass from the center of the density peak 
for several LTC models. The enclosed mass is defined as 
the mass within a sphere of a given radius. The mass 
accretion rate M is defined as the average of M over 
particles: 



M 



1 



™shcll 



i=l 



4irpi (vi - Dp k ) • (x pk - a;i)|a;pk - xi 



. ( 3 ) 

where n s h e ii denotes the number of SPH particles within 
a thin shell (r, r + <5r), and the summation is taken over 
the particles within the shell. a; p k and t> p k denote the 
position and velocity of the particle at the density peak, 
and Xi and Vi are those of the particle in the shell. 

Five models in which the core collapses successfully 
are plotted. Four models are LTC models with ionizing 
photons, while a model denoted by a green curve rep- 
resents a UV-free (no feedback) case as a reference. In 
the case of no feedback, the accretion timescale intersects 
the KH contraction timescale for the enclosed mass of ~ 
6 x 10 2 M(7). This is comparable to the previous estimates 
dBromm. Coppi fe Larson|[l 999: Na kamura fc Umemural 
120011 lAbel. Bryan. Normanll2002t lYoshidall2006t )~ 

On the other hand, with radiative feedback, the accre- 
tion timescale becomes longer than that in no feedback 
case. For instance, in the model with n on = 3 x 10 3 cm -3 
and D = lOOpc (blue dotted curve), the mass accreted in 
the KH timescale is reduced to a level of ~ 20A/q. There- 
fore, the radiative feedback has a significant impact on 
the mass of secondary Pop III stars. We also remark 
that this effect has strong dependence on n on and D. 
Blue and red dotted curves in FigfTTI denote the model 
close to the boundary of successful collapse region (see 
Fig. [9]), whereas blue and red solid curves correspond to 
the models slightly distant from the boundary. The mass 
accretion timescale in latter two models is also prolonged 
by radiative feedback, but the effect is not so dramatic 
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Fig. 4. — Hydrogen molecule distributions on a plane including the symmetry axis are shown at lMyr and 2Myr after the UV irradiation 
for HTC-ION and LTC-ION models. The source star is located on the left boundary of each panel. H2 fractions are shown by colored dots 
according to the color legend. 

as the former two models. In these cases, the expected 
mass of secondary stars is as large as « WOMq. 

7. CONCLUSIONS AND DISCUSSION 

We have explored the evolution of a collapsing cloud 
nearby a hrst luminous star to examine quantitatively ra- 
diation hydrodynamic feedback on secondary Pop III star 
formation. To investigate this issue, three-dimensional 
radiation SPH simulations have been performed, where 
the radiative transfer of ionizing photons and H2 disso- 
ciating photons from a first star is self-consistently in- 
corporated. We have derived the numerical criteria for 
the subsequent Pop III star formation. As a result, we 
have found that an H2 shell formed ahead of an ionizing 
front works effectively to shield the H2 dissociating ra- 
diation, leading to the positive feedback. On the other 
hand, a shock associated with an ionizing front has a pos- 
itive effect and a negative effect. The positive one is the 
compression and heating of gas cloud that promote H2 
formation and thereby accelerate the gravitational insta- 



bility of the cloud. The negative one is the stripping of 
gas envelope that reduces the mass accretion rate, lead- 
ing to the formation of less massive Pop III stars. In 
a small core, the compression and heating by a shock 
turns to be devastating to hinder the cloud collapse. In 
comparison of the mass accretion timescale with Kelvin- 
Helmholtz timescale, we have found that the mass of 
secondary Pop III stars could be reduced to rj 20M©. 

The possibility that the formation of low-mass Pop III 
stars with < 100M Q is driven by radiation-hydrodynamic 
feedback would be important in the environment where 
multiple stars can form within a ~ (lOOpc) 3 volume. 
In the conventional one-star-per-halo picture, this low- 
mass star formation mode would not be so impor- 
tant. However, in the case of first galaxies, this 
mode could play an important r ole as well as the 
other modes driven by HP cooling (TUehara fe Inutsuka 
2000 1; iNakamur a fc Umemur a 2002; Naga kura fe O mukai 
200a I Johnson fc Bromml 120061 : iGreif fc Bromml 120061 : 



Fig. 5. — Same as Fi 
depicted in units of cm - 




LCT-ION 2Myi 



10^ 



HCT-ION 2Myr 



except that the colors of the dots represent the number density of hydrogen nucleous. The color legend is 



lYoshida. Omukai fe HernquistJ 12007). In addition, even 
in a minihalo, the assumption of one-star-per-halo it- 
self could be modified by resolving dark matter evolution 
with very high accuracy (|Umemura et al.ll2009f L 

If such lower-mass Pop III stars are formed, the 
radiation-hydrodynamic feedback on further star forma- 
tion could be changed, since the importance of ionizing 
radiation relative to H2 dissocia ting radiation is differ- 
ent from that of a 120M Q star. lHasegawa et al.1 (|2009h 
have investigated the dependence on the stellar mass 
of radiation-hydrodynamic feedback and concluded that 
the ionizing radiation from ~ 25Mq stars does not have 
significant impact on the neighboring star formation, but 
the feedback is essentially determined by H2 dissociating 
radiation, which causes only a negative effect. 

In the present simulations, we omitted a few processes. 
First, we do not include the HD chemistry/cooing. It is 
well known that HD is an important coolant in case once 
the gas is heated above 10 4 K. In that case, H 2 cooling 
proceeds faster than the recombination of hydrogen atom 



at several thousand Kelvin, leadi ng to a large H2 fraction , 
because of "frozen" electrons (IShapiro fc Kane 1 119871 : 
iSusa et all fl998t lOh fe Haiman I 12002ft. HD molecules 
are formed from these abund ant H2 molecules as the 
collapse of cloud proceedsfe.g , lUehara fe Inutsukal 120001 ; 
iNakamura fc Umemural l2002f h Thus, HD cooling could 
be important when the ionized gas recollapses after the 
death of the source star. However, in the present sim- 
ulations, the collapsing density peaks do not experience 
such a hot phase (FigJ2j) , because of the self-shielding of 
the core. The shock associated to the I-front is also too 
mild to ionize the material (again FigfS]). Thus, as far as 
the collapse criteria of the core is concerned, HD cooling 
do not play an important role. Secondly, we also do not 
includ e the H2-H + co ll ision- induced H2 cooling. Accord- 
ing to iGlover fc Abell (|2008ft , H2-H 4 " collision dominates 
the excitation of H2 in case electron fraction y c > 10~ 2 
at T ~ 10 4 K, and y c > 10~ 3 at T ~ 10 3 K. These con- 
ditions are marginally satisfied when the gas cools down 
from ionized state. Thus, H2-H+ collision-induced H2 
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Fig. 6. — Same as FigU except that the colors of the dots represent the temperature of hydrogen nucleous. The color legend is depicted 
in units of Kelvin. 



cooling can play an important role, when the ionized gas 
recollapses after the source star turned off. However, 
again in our simulations, it is not so important, since the 
gas at the density peak is not ionized in its history. It 
is also worth noting that cooling and dynamics in the 
neighbo r of H2 shell could be aff ected by these cooling 
process ()Whalen fc Norm an 2008), however, it is unlikely 
that this effect changes the present results. 

We are grateful to all the collaborators in Cosmolog- 



ical Radiative Transfer Codes Comparison Project for 
fruitful discussions during the three workshops. Numeri- 
cal simulations have been performed with computational 
facilities at Center for Computational Sciences in Uni- 
versity of Tsukuba and with computational facilities in 
Rikkyo University. This work was supported in part by 
the FIRST project based on Grants-in-Aid for Specially 
Promoted Research by MEXT (16002003) and Grant-in- 
Aid for Scientific Research (S) by JSPS (20224002) and 
Inamori Research Foundation. 
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